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ON : 

0\ ■ Abstract 

; 

. We analyse the Structure Function collaboration data on the quark pseudoscalar vertex 

' and extract the Goldstone boson pole contribution, in l/p^. The strength of the pole is 

^ . found to be quite large at presently accessible scales. We draw the important consequences 

'"^ ' of this finding for the various definitions of quark masses (short distance and Gcorgi-Politzer) , 

\l , and point out problems with the operator product expansion and with the non-perturbative 

' renormalisation method. 

(N ■ 
> ■ 

QQ I 1. Continuum model for the quark pseudoscalar vertex 

vn ■ 

■ It is well known that the quark pseudoscalar (PS) vertex contains a non-perturbative contri- 

, bution from the Goldstone boson, in the continuum. On the lattice, the use of a non perturbative 

renormalisation scheme |^] makes this contribution manifest, although it should go to zero for 
Q^ ' large momentum transfers. The purpose of this letter is to extract it from lattice simulation 

data, and to show that it is not negligible for presently accessible scales. In particular, it must 
be subtracted when evaluating the short distance quark masses from the lattice through the 
pi ' non perturbative method of ref. This method uses the off-shell axial Ward identity (AWl) 

^ . and renormalises the mass in the momentum subtraction (MOM) scheme of ref. Q , which can 

afterwards be related to the MS scheme. The MOM renormalisation involves the pseudoscalar 
vertex, whence the necessity of the subtraction of the Goldstone contribution to extract the 
^ ■ short distance quantity^. For physical u, d quarks, the Goldstone contribution becomes very 

large, larger than the perturbative part; this corresponds to a very large dynamical u, d mass, 
larger than the usual current mass at the scales accessible to standard lattice calculations. 

The expected behaviour of the pseudoscalar vertex in the continuum has been described as 
follows in the 70's in the works of Lane and Pagels ^ and others. Near the chiral limit, the 
one-particle-irreducible PS quark vertex A5 can be described through a perturbative contribution 
plus a non-perturbative Goldstone boson contribution. 

The perturbative contribution is of course 1 x 75, with QCD radiative corrections, which 
according to the renormalisation group lead to a logarithmic behaviour [^^(p^)]^/^^ for Np = 0. 

As to the non-perturbative contribution, firstly, according to PCAC, the Goldstone boson 
must dominate other pole contributions in the PS vertex near the chiral limit: the coupling 
of the pion to the PS vertex indeed gives a pole ~ ^/{o^ + w-^); where q is the momentum 
transferred at the vertex; other poles (radial excitations) contributions are suppressed in the 
chiral limit. At g = 0, in terms of the quark mass m, this gives a 1/m pole contribution. We 
emphasize that this pole contribution explodes at = in the chiral limit, i.e. it is singular 
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^Laboratoire associe au CNRS-URA D00063. 

^ The problem would be quite different if other methods are used to extract quark masses, see e.g. ref. 
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on the border of the "physical" region of momenta. This well-known fact seems to have been 
underestimated, especially in its consequence for the renormalisation of mass on the lattice near 
the chiral limit (see below, section 3). 

Secondly, according to the Wilson operator product expansion (OPE), the non-perturbative 
contribution must be power behaved, i.e. at large it must drop as l/p^ up to logs. This power 
behavior, as shown by Lane, is dictated by one of the pion Bethe-Salpeter amplitude in the OPE. 
In the latter, the dominant operator is the vacuum to pion matrix element of the pseudoscalar 
density. This leads to by the canonical dimensions and further to logs, [as(p^)]'^/^^ at 

Finally, a very important relation, emphasized by Lane and Pagels and derived from the 
Ward identity, connects directly the forward pseudoscalar vertex to the scalar part of the prop- 
agator (they are essentially proportional)^. Hence the study of propagator OPE by Politzer |^] 
closely parallels the above considerations on the PS vertex, with the dominant non-perturbative 
(power) contribution corresponding to the quark condensate, which corresponds to the Gold- 
stone pole contribution in the PS vertex. We postpone the discussion of the propagator lattice 
data, which involves delicate problems of improvement, to another paper [^]. Nevertheless, in 
the last part of the letter, using only the vertex lattice data, we draw important consequences 
for the propagator thanks to this Ward identity. 

In addition, we use the Ward identity to borrow the two-loop corrections to the above 
logarithmic factors, from those calculated in the case of the propagator by the Rome group 
and by Pascual and de Rafael |^]. 

In the following, we will show that the lattice data behave as expected from the continuum 
theoretical expectations. To relate lattice to physical numbers, we will use the lattice unit 

= 1.9 GeV at /? = 6.0 and aAj^g = 0.174 for the QCD scale parameter in the quenched 
approximation. We must stress that nothing in our discussion depends critically on the precise 
values of these parameters: this work is not oriented towards accurate numerical determinations, 
but rather towards questions of principle. 

2. Fit of lattice data for Fs as function of the hopping parameter k 

Recently |10|, Paul Rakow presented, on the behalf of the Structure Function (QCDSF) 
collaboration, the bare vertex F5 data at /? = 6.0 (Fig. 1 of ||lO[| ) through the product 

Iq f 5(p^) = C X amq 



amqf5{p^) = Cx amq F5(g = 0,p^) (1) 



as a function of a^p^ at three k values, with: 

^ _ Tr (75 As) 



2) 



4 

This is the only contribution that survives at g = in the continuum limit. By construction, C 
is a constant with the value [10|: 



C = 0.75, (3) 

■'in the work of Lane, as well as in some subsequent works, one quotes another power of logs: [as(p^)]^*''^^. 
This latter power corresponds to the anomalous dimension of the local operator tp^stp. But one requires a gluon 
exchange at lowest order, which gives an additional factor. This was shown by Politzer Q for the corresponding 
^ptf) contribution to the propagator 

''Let us also recall that in the context of the lattice, this relation has been discussed and used by the Rome 
group to fix Za on quark states M. 
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Figure 1: (a) The value of the coefficient of a^p'^A{p'^) in eq. @ from the lattice data extrapolated 
at Kc, (b) Our fit to the extrapolated data. 



defined in |10] so that the numbers in r.h.s. of eq. ([l|) equate approximately the lattice data for 
the scalar part of the bare improved propagator (same Fig. 1 of ||l^) on some range of momenta, 
increasing with k. This approximate equality corresponds to the Ward identity to be discussed 
in section 4, but we are not concerned with it, as we rely only on the vertex data. The QCDSF 
data were obtained with the SW improved fermionic action [|ll| at csw = 1.769, Np = 0, in the 
Landau gauge, and p"^ denotes the following lattice definition of the momentum squared. 
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where a is the lattice unit. This notation will be used herefrom, and identified with the continuum 
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p". We take a 
the bare mass 



Let us write 



1.9 ± 0.1 GeV at /? = 6.0 |12]. We also use the standard lattice definition of 



am„ 
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amgf 5(p^) = A{p^) + amqB{p'^) 



(5) 



(6) 



We shall first perform an extrapolation linear in ruq of the three datasets to the chiral limit at 
lie = 0.1352, for each value of p^. This fit gives us both A{p'^) and B{p'^), as shown in Figs. 1 
and 2. 

As T^iq = 0,p2) cx [A{p^)/amg + Bip"^)], it seems reasonable to identify the A{p'^) term as 
the Goldstone contribution (i.e. as the pole in aniq), and the second one as the perturbative 
contribution, if we are sufficiently close to the chiral limit. The continuum Ward identity enables 
us to borrow the two-loop renormalisation group improved factors, from those calculated in the 
case of the propagator by the Rome group ^ for the perturbative part: 



B{p') = Bo X [a,ip')] 



2^14/11 



1 + 8.1 



Qs(p^) 
47r 



(7) 
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Figure 2: (a) The value of B{p'^) in eq. @ from our extrapolation at Kc, compared with 
6q- (0); (b) Our fit to the lattice data in the region of high for Bq = 1.735. The plain curves 
correspond from top to bottom respectively to aruq = 0.148, 0.078, and 0.028, and the dashed 
one to the Goldstone boson contribution. 



and by Pascual and de Rafael 



for the non-perturbative part: 



An X 



2^l7/ll 



1 + 22.0 



4tt 



(8) 



with Aq and Bq some constants. These two-loop corrections are valid for the MOM renormalisa- 
tion scheme, in the Landau gauge, but with Ug taken as the MS coupling constant^. As already 
mentioned, we take aAj^g = 0.174 at Np = from the three gluon coupling measurement with 
asymmetric momenta 14]. Note that, a priori, the evolution formulae properly apply to the 
renormalised propagator, which is the bare one divided by Z^(p^) 0; but, at least theoretically 
and in the continuum, Z^{p'^) should evolve very slowly, since 7^ = at one loop in the Landau 
gauge, and the two-loop correction seems also to have very little effect Lattice data [||, ||] 
confirm this perturbative argument. Then, we conclude that the two-loop corrections should be 
obeyed by the bare PS vertex with good accuracy. 

The lattice data turn out to be quite close to the continuum theoretical expectations near 
the chiral limit, confirming the above interpretation of A and B. Indeed one finds the following: 
• A{p'^) is behaving remarkably close to 1/p'^ over a large interval of p^, see Fig. 1(a). From the 
numerical analysis, the Goldstone contribution appears to be very large. Indeed, 



aV^(p^) ^ 0.015 



(9) 



from the lowest point a^p^ = 0.16. 

On the other hand, we do not see the log factors expected from the perturbative calculation, 
eq.(^), which remains to be understood; we shall discuss further problems in the perturbative 
evaluation of the Wilson coefficient at the end of the paper. 



^for which we use the expression 47r/as(g^) = 11 log 



In this letter, we use the convention of ref.l 



MS 

for the Z's. 



^ log log 
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• Bijp') is found to evolve closely to ln(p^/AQ^£,) ^/^^^ more precisely it is evolving in good 
conformity with the two-loop MOM renormalisation formula quoted above. We obtain 

Bq = 1.735 (10) 

which provides a very good fit to the data for larger than 2 GeV^, see Fig. 2(b). 

The Goldstone contribution is felt already at rather large quark masses and, for physical 
It, d quarks it is in fact very large: A{p^) is larger than the perturbative part amqBljP') even at 
rather large p^, as shown in Fig. 3(a) and further discussed in section 4. 

3. Consequences on the short-distance mass from the lattice 

Let us then recall that the pseudoscalar vertex is an important ingredient in the method 
first developed by the Rome group to determine the short distance quark masses. One starts 
from the "axial" MOM renormalised quark mass given by: 



am^^^f "(/.2) = p^A_ (11) 



where p \s a. dimensionless parameter determined from a ratio of matrix elements involving the 
pion and the bare axial and pseudoscalar bilinear operators. Za is the standard renormalisation 
of the axial current, determined from the axial Ward identity or approximated through the MOM 
renormalisation non-perturbative method, or from one- loop perturbation theory. Finally Zp{fi'^) 
is defined through the MOM renormalisation condition for the vertex and is determined by the 
same non-perturbative method from the PS quark vertex at q = 0, or again approximated from 
one-loop perturbation theory]^ Then, from m^^j""(^^), one can deduce the standard short 
distance masses, for instance the MS mass, but this requires to work in the perturbative, short 
distance, region. 

Then our findings for the PS vertex have important consequences. Let us note that, in 
principle, the fact that at g = 0, the PS vertex not only is influenced by a large Goldstone 
contribution (as noted in |^], |]l8|), but really explodes in the chiral limit, is crucial for the 
MOM procedure of renormalisation, since the renormalisation constant Zp is defined as: 

ZMOM^^2^ = (12) 

This definition ensures that 

rf(g = 0,/) = 1 atp2 = ^2 (-^3) 

and one concludes from it that Zp has a trivial chiral limit at fixed p^. Indeed, near Kc, since 
cx nig and T^iq = 0,p^) oc l/rUg, tends to zero, or its inverse goes to infinity. 

3.1. Calculation of Zp 

One can translate the above fit of Fs into an expression for Zp'-'^ , or rather its inverse 
which is more directly physical: 

1 T^ip") Azip"^ 



+ Bz{p') (14) 



^Of course, using the PS quark vertex for normalisation is a particular way to define the MOM scheme, which, 
in the non-perturbative regime, is not necessarily compatible with e.g. the use of the scalar vertex. See remarks 
below. 
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Figure 3: (a) The value of anigT^^q = 0,p^) for light quarks; (b) The value of the dynamical 
u, d masses. 



where ^^(p^) = Aip") /[C Z^^p")] and Bzip") = B{p^)/[CZ^{p% 

This is to be contrasted with usual fits, which assume that Zp is linear in aniq of eq. 
One notes that since Z^ljp') is weakly dependent on and on k, and close to 1, the expression 
is quite similar to the preceding one: it consists in a first term which is approximately in l/ruq, 
and the second one which is approximately constant^. Just as for the Ts vertex, the former 
corresponds to the non perturbative Goldstone contribution while the latter, Bz, corresponds 
to the short distance contribution, as measured by the lattice numerical simulation, including 
all the orders of perturbation theory by a non perturbative method. To give numbers, we need 



now values for Z^, which we borrow from the data of the Rome group |16|, with an improvement 
procedure trying to parallel as much as possible the one followed by Rakow for the scalar part. 
However, we would like to emphasize that all the qualitative conclusions of this paper are 
independent of precise values of long as they stay around 1. 

At a'^p'^ = 1 which is close to the standard reference point p = 2 GeV, and at k = 0.1342, 
which is the k closest to the chiral limit, we estimate: 



Z^ ~ 0.85 



One has then numerically: 



Az ~ 0.023, 



Bz ^ 1. 



(15) 



(16) 



We emphasize that this identification of the two contributions does not rely on Boosted Per- 
turbation Theory (BPT), but purely on numerical simulations. We shall refer to Bz as the 
"short distance contribution" not to be confused with its one-loop perturbative estimate. Now, 
we shall compare to BPT estimates here and in the following for illustrative purposes only. 



*D. Becirevic has now done a fit along the above lines with the data of the Orsay-Rome group and found 
roughly similar conclusions. 

^One could be worried by the fact that the dependence of Z^p on k could generate from Az/[{1/k. — 1/kc)/2] 
an additional contribution to Zp^ independent of k. However, it can be seen that since Az decreases rapidly, this 
contribution is not large with respect to the one coming from Bz- 
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Bz is indeed very close to the one-loop standard BPT evaluation Zp^{a?p'^ = 1) = 1.7 (from 
Zp = 0.59) in the chiral limit, with = 1.6822, whereas the non-perturbative term contributes 
around 0.023/amg ~ 0.8 on a total of 2.7. Therefore, the departure of Zp from its one-loop 
perturbative evaluation seems to be essentially due to the Goldstone boson contribution. Higher- 
order radiative corrections do not seem to be very largeP^. However, these conclusions could be 
sensitive to details of the data or of the fits. 



3.2. Calculation of the MS mass 

We can now convert our results for the renormalised mass m^^^"" into a calculation of the 
MS mass. In ref. [^], it is suggested that the use of the MOM non-perturbative determination 
of Zp (with linear extrapolation in n) in the AWI method improves the results for the MS mass 
with respect to previous determinations of Zp by one-loop perturbative calculations. 

However, in principle, to make the conversion to a short distance mass, we must still make 
sure that we work in the perturbative regime. Now, we have isolated the Goldstone boson 
contribution which is essentially non-perturbative as it does not correspond to higher order 
contributions but rather to power corrections. The fact that the non-perturbative estimate of 
the full Zp differs sizeably from the short distance Bz already at the measured kappas, is a 
signal that it is not presently possible to work at high enough for the Goldstone contribution 
to be negligible. Hence, we must first subtract it from Zt} . 



ZSuUr^pZ^ -= Zp^(p^) -'-l±^ = Bz{p') (17) 



Az{p 



arun 



Numerically, the remaining short distance Bz is in fact close to the one-loop BPT estimate of 
Zp, as we have just seen. Indeed, at a^p^ = 1, with this subtraction and using again Z^ = 0.85 
near the chiral limit, we find Zl"**'' = 1/1.88 = 0.53 which corresponds to the fully resummed 
short-distance contribution determined directly from the lattice dataQ To get an estimate of 
the consequences on the light quark masses, we use [^2| (which uses the notation am for p) 



where one finds p ~ aruq, with aniu^d = 0.001836, and we take Za = 0.79 [23|. We then find: 

a m^7'^''" ~ 0.0027 (18) 
Converted into the MS scheme through pi] 



16 a, 



^M.^^Lanaau (^l--^j (19) 

with asia^p^ = 1) = 0.25 at two loops, this gives: 

a m^J ~ 0.0024 (20) 

therefore about 4.6 MeV at Np = 0. One would obtain about 6 MeV if one used the full MOM 
non-perturbative estimate of Zp linearly extrapolated to the chiral limit and 4.2 MeV from one- 
loop BPT. Note that these numbers are only indicative; in view of the many uncertainties in 
the subtraction procedure, we do not try to discuss the other sources of error necessary to give 



^°These two possible explanations were suggested in ref. Q. 

We take Zp"''''" approximatively independent of mass, except for a small variation of Z^. One would otherwise 
require a fit of m^Fs with one more term in m,. 



7 



a real determination of the mass. Our aim is only to underline the necessity of the subtraction 
of the Goldstone contribution. 

Despite the fact that in this case the two methods lead to comparable results, the subtraction 
method just described, which does not rely on perturbation theory but which rather uses the 
lattice data directly, is in general superior to that based on the BPT estimate, and knowledge 
of the full Zp'^^ is in general necessary even if we aim at measuring short distance quantities. 

Indeed, the BPT method has the following drawbacks. Firstly, the unknown higher-order 
perturbative corrections may be large. Secondly, the one-loop perturbative evaluation can be 
tadpole-improved in many ways, potentially leading to very different estimates. Finally, the one- 
loop estimate, even if tadpole-improved, does not automatically follow the behavior dictated by 
the renormalisation group; the problem is then cured by taking the one- loop estimate at some 
momentum, and then imposing the renormalisation group evolution for the other momenta, but 
this is obviously presenting a rather arbitrary choice of a privileged point. 

On the other hand, the subtraction method used here avoids these problems. It amounts in 
general to a non-perturbative measurement of the Z's, followed by the evaluation and removal 
of the pole contributions. As shown above in Fig. 2(b), this procedure leads to a result which 
evolves per se according to the renormalisation group. 

The evaluation of the pole contribution from the pion is especially easy, because of its 
particular singular nature at ruq = 0. Other pole contributions are expected to be smaller, 
because they are regular in the chiral limit. If one were to subtract them, one could only rely 
on the expected power behavior of the particle vertex function, and their extraction would thus 
be more difficult. 

4. Consequences for the renormalised mass a la Georgi-Politzer 

In this section, we show that the Goldstone boson contribution to the PS vertex, which 
is only parasitical in the calculation of MS masses, and has to be subtracted as shown in the 
previous sections, retains an important physical meaning, as can be seen through the use of 
other definitions of renormalised quark masses. 

4.1. Physical relevance of the full renormalised axial mass 

One should remember that the full axial MOM renormalised mass, which is calculated 
through: 

am^A^t^ip^) = pZA[Zpip^)]~' (21) 

with ZpijP') not submitted to the above subtraction, retains a physical significance by itself, since 
it is the mass defined through the divergence of the axial current, with the corresponding natural 
renormalisation condition which consists in setting the pseudoscalar vertex to 1 on quark 
states at the renormalisation scaleP). We stress that in contrast to the standard MS current 
mass, it does not vanish in the chiral limit, because the chiral limit p ^ is compensated by the 
pole in Zp . The meaning of this last result will be further developed in the next subsection. 
Of course, it is unpleasant to have renormalised quantities with a rather queer behavior in the 

^^Recall then that the scalar vertex should be consistently normalised through the Ward identities, as repeat- 
edly emphasized by the Rome group, therefrom one sees that Zs/Zp must be independent of p^, in obvious 
contradiction with imposing simultaneously a similar MOM condition for the scalar vertex. 
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chiral limit. Indeed, for instance, hadronic matrix elements of the renormalised pseudoscalar 
density which do not have a pion pole should tend to zero in the chiral limit- and therefore also, 
the ones of the scalar density defined in accordance with the Ward identity. But as we hope to 
have shown, this is an unavoidable consequence of the non-perturbative method as applied to 
the pseudoscalar density. Of course, one could avoid this by preferring the corresponding MOM 
normalisation condition for the scalar vertex. But then, one must remember another aspect of 
the axial MOM renormalised mass, which gives it another important physical significance, and 
which we shall now discuss. 



4.2. Relation with the mass as defined by the propagator 

It can indeed be easily shown that the "axial" renormalised mass is also essentially identical 
to a standard renormalised mass defined through the scalar part of the propagator, as becomes 
obvious through the renormalised axial Ward identity at zero transfer (this has been recalled in 
the talk of Rakow); in fact it is essentially identical to the Georgi-Politzer mass. Let us indeed 
write the Ward identity]^: 

mi^^f«(^^)rf (, = 0,/,^2) = (22) 

Here, it is assumed of course that the propagator is also calculated in the Landau gauge, and 
that the renormalisation of the vertex and of the propagator are both consistently performed 
according to the MOM scheme. Therefore the propagator is normalised in the Euclidean region 
according to^: 

S],'{p,fi^)=i ^ + m'l''{fi^) at /^m' (23) 
which is nothing else than the Georgi-Politzer renormalisation condition. 
Setting = /u^, one obtains from eqs. (p^, ( p2D and ([2^): 

^AWI (M j - Ip2=;,2 - [fi ) (24) 

therefore the announced identity of m^p|)j""(/i^) and the Georgi-Politzer mass function m^^^j?) 
is derived. 

Through eq. (p^, the l/p^ power contribution in Zp^ corresponding to the Goldstone 
boson is related to a similar contribution in the scalar part of the propagator, which represents 
a dynamically generated mass for light quarks, though off-shell, gauge-dependent and Euclidean. 



This is a well-known |24| signal of the spontaneous breakdown of the chiral symmetry. We can 



then translate our knowledge of the PS vertex into information on this dynamical mass. 
4.3. Physical consequences 

Numerically, for the non-perturbative contribution, which corresponds to the chiral limit of 
m^^f"", one finds at a'^p^ = 1, from eq. (21), with, as before, Zp from eq. (0), Za = 0.79 and 



^•^Here we use eq. ( |22[ ) as a constraint allowing to calculate the r.h.s. from the l.h.s. We do not require the 
independent input of the scalar part of the propagator data. As a side remark, as a check of the validity of eq. ( p2[ ) 
on the lattice QCDSF propagator data, we note that it would imply for the C defined above C = ZAp/amq, which 
is not far from their C — 0.75. 

We disregard here the difference between this standard MOM condition for and the derivative MOM 
condition derived from the vector Ward identity by the Rome group |^], which seems to lead numerically to very 
small differences. 
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a m%^{a^p^ = 1) ~ 0.018 (25) 

therefore around 34 MeV at p = 1.9 GeV. 

We have obtained analogous resuhs directly from the scalar part of the bare improved 
propagator as given in [|^. Extended considerations and estimates on the propagator wih be 
given in a forthcoming paper [^. 

At large p, we want to stress that the non-perturbative contribution to the mass remains of 
the order of the u,d perturbative masses and even larger, see Fig. 3(a). It must be emphasized 
that this result is rather safe, since it does not depend critically on improvement procedures. 

The sign is as expected, but the magnitude is much larger than what one would expect from 
the estimate by the qq condensate and a perturbative calculation of the Wilson coefficient [^: 

GPr 2. 4 < 0|gg|0 > 

mR [p ) = --vTQs 2 (26) 

3 

Indeed, with a^^{a^p^ = 1) ~ 0.3 (one-loop), and with a standard MS renormalised evaluation 
of the qq condensate of —(225 MeV)'^, at = 1 GeV , rescaled at aP'p^ = 1 by a factor 1.18, one 
would find an answer lower by a factor ten. However, a large value is inevitable if we admit as 
usual that at moderately low the mass is of the order of a constituent mass, i.e. several tens 
of MeV, and if we take into account that the decrease is only in 1/p^, as found with remarkable 
accuracy on the lattice data of the QCDSF group. 

That a value of several tens of MeV is needed at low p^ is confirmed by direct lattice 
calculations of the propagator, to be compared with our estimate of Fig. 3(b): 
• in configuration space, the propagator in time S{t) has been found exponential over a large 



range of t, with a coefficient of the exponential around 300 MeV |17, Pq ]; 

• in momentum space, at the lowest points, corresponding to p^ — > 0, where it should be 
insensitive to the improvement, the scalar part of the inverse propagator is found to be around 



similar values of 200 — 400 MeV [10, for the lowest mass k = 0.1342. A large value for the 
coefficient of 1/p^ was also found in a phenomenology of the pion as a Goldstone boson [p^ ], 
based on an assumption: 

„gP . = (27) 

where m£, is a free parameter. The phenomenology seems to require m£, ~ 300 MeV, m^^{p'^) ~ 
27 MeV aX p = 2 GeV, in agreement with the above estimate. 

It must be emphasized that this large non-perturbative contribution does not contradict di- 
rectly the sum rule calculation of correlators, which uses a perturbative evaluation of the quark 
propagator. Indeed, as in the calculation of propagator by Politzer, in sum rule calculations 
non-perturbative contributions are consistently added through condensates. Moreover, as em- 
phasized by Pascual and de Rafael |^, the condensate contributions are quite different for the 
quark propagator and the correlators, therefore our finding for the full quark propagator does 
not have a direct impact on sum rules. 

Of course, the discrepancy with the naive perturbative estimate of the Wilson coefficient is 
nevertheless worrying and deserves further reflexion since the latter is known to work in general. 
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One could imagine that the perturbative calculation of the Wilson coefficient is not valid for 
some particular reason. In this direction, one must observe that the two-loop correction is very 
large, around 50% of the first order at a?p'^ = 1. It may be therefore that the perturbative 
expansion happens to fail. 

Conclusion 

The lattice numerical calculations can be seen as the triumph of the general theoretical pre- 
dictions of the 70's for the quark pseudoscalar vertex. Nevertheless, a striking and unexpected 
feature of the lattice data is the very large size of the Goldstone boson contribution to the pseu- 
doscalar vertex. This corresponds, through the Ward identity, to a very large non-perturbative 
contribution to the renormalised mass function of Georgi and Politzer, by far larger than what 
is expected from the quark condensate and a perturbative evaluation of the Wilson coefficient, 
but in agreement with other physical expectations. 

These large non-perturbative contributions then give a warning that there may be a possible 
problem with the use of lowest order perturbation theory in the estimate of the Wilson coefficient 
of condensates. This question requires more investigation. 

The Goldstone contribution must be subtracted from the pseudoscalar vertex to calculate 
the short distance mass from a normalisation of this pseudoscalar vertex. This has important 
numerical consequences. The short distance Zp to be used should be sizeably larger than the 
one measured directly on the lattice. 

Finally, it must not be forgotten that lattice artefacts may still be large at /3 = 6, , and that 
the vertex has not been improved with rotations as an off-shell Green function, therefore one 
can expect large uncertainties on the quantitative estimate of the effect at large p^. 

Note added in proof: When completing this paper, we became aware of the work of the JLQCD 
collaboration, presented at Denver conference, where certain parallel conclusions on Zp have 
been drawn from lattice staggered fermion data. Paul Rakow has also drawn our attention to 
ref. Jl^/ , where connected observations on PS vertex and the scheme of Pagels are made. 
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